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ABSTRACT 


The principle of linear superposition is applied to Bleich's 
two general solutions pertaining to hull-deckhouse interaction. The 
meometing longitudinal stress distributions for three separate mathemati- 
cal models are compared with stress distribution results obtained 
through the use of finite element techniques. In reference to Bleich's 
theory. procedures are derived for the determination of the stress dis- 
tribution at locations away from the center of the hull-deckhouse struc- 
ture. Anaiytic and finite element procedures are also described for 
the determination of the deck stiffness. 


The principle of linear superposition is found to be valid 
for Bleich's theory. Because the theory neglects shear lag effects, 
its application is recommended only for the center portion of struc- 
tures with relatively long deckhouses. 
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INTRODUCTION 


In 1953 H. H. Bleich presented a paper in the "Journal of Applied 
Mechanics" entitled "Non-Linear Distribution of Bending Stresses Due 
to Distortion of the Cross Section." In this paper he derived a via- 
ble analytical solution to the problem of hull-deckhouse interaction. 
However, because of the difficulty involved in evaluating the two- 
coupled differential equations for real life situations Bleich was 
only able to present two simplified solutions. One solution to the 
differential equations was simplified by considering only bending 
moment loads and neglecting the distributed loads. The other solution 
considered only distributed loads and neglected bending moment loads. 
Although these simplifying assumptions are not justified in any real 
Ship, the methods of solution are, in actuality, straightforward and 
rather simple to obtain. The proposal is that the principal of linear 
Superposition can be applied to Bleich's solution in order to deter- 
mine the total solution. Even though the title of Bleich's article 
implies non-linearity, the calculated stress distributions found were 
actually linear in form with the only discontinuity or break annpearing 
at the main deck level. In order to confirm or disprove this hypothe- 
sis, finite element methods using the ICES STRUDL program developed 
at M.1I.T. will be employed to find the total bending stress distribu- 
tion. A comparison of the linear superposition of the two solutions 
using Bleich's method and that of the total ICES STRUDL solution should 
give some indication of the applicability of linear superposition to 


Bleich's method. If superposition if applicable, a rather simple 
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analytical solution could then be available for the determination of 


the actual total bending stress distribution. 





CHAPTER I ~- BLEICH'S TWO-BEAM THEORY 


Consider the problem of two separate beams forced to act together 
by horizontal shear forces and vertical forces acting at the junction 
of hull and deckhouse. In essence, the system consists of a beam 
elastically supported by another beam with a shear connection to enforce 
equal strains at the deck level, Figure 1.1 and 1.2. The vertical 
forces are due to elastic resistance of the deck framing or bulkheads 
against the motion of the deckhouse with respect to the hull. Navier's 
hypothesis for a structure acting as a beam is applicable only if the 
cross sections of the beam do not distort. Since the side plating of 
a deckhouse is not coplanar with cag ships side,the vertical deflections 
of the two are not necessarily the same. Navier's hypothesis that 
plane sections remain plane and the conventional theory of bending of 
beams are not applicable for the entire structure. For this analysis, 
however, Bleich assumed that the symmetrical cross section consists of 
two non-distorting portions which can move relative to each other and 
that Navier's hypothesis is applicable for the hull and deckhouse each 
by itself. 

In the structure shown in Figure 1.1, the lower hollow box beam 
represents the hull and is of length L while the upper box, the deck- 
house, is shorter and is of length 1. Both boxes are assumed to be of 
constant cross section. The cross-sectional area and the moment of 
inertia of the deckhouse and of the hull are Al. I: and Ay. I,. 
respectively, and the distances of the respective centers of gravity 


from each other and from the deck are a, O42, anc 15a. 
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mememplified Analysis of Twortell Structure 

In this analysis the assumption is made that the deck A-B, kigure 
1.1, and its supnvorts have no stiffness and, hence, will not resist 
any relative vertical movements between hull and deckhouse. This 
simplification is not justified for any real ship system and is dis- 
carded in the final analysis; however, some important results are 
derived. 

At a distance z in the free body diagram of Figure 1.2, the moment 


and direct forces in the deckhouse and hull are Mi Nae and My Nos 


respectively, with positive moments producing compression at the top 


of the respective units. Direct forces Ny and N, are positive if they 


create tension. The external loads acting to the left of the section 
have a moment M. A shear force T of unknown magnitude will act on the 
underside of the deckhouse, and a similar force T will act in the oppo- 
site direction on the hull. Equilibrium requires that: 

Ny 1 7 7 ae,T (lia) 


N, = T, M, =M- aot ° (1b) 


th 
t 
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Since it is assumed that Navier's hypothesis is valid for deckhouse 
and hull separately, the longitudinal stress distribution across a 


section can be determined for the deckhouse 


et. ts. (2a) 
oe 


andecor the hull 
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9 (2b) 


with tension stresses counted as positive. 





At the junction of the house and hull, the stresses O4 and C5 


must be equal, with x) = 7 ad) and Xo = ado. Using Equations 1 and 


Peeeene value T is found to be: 
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1 | 
T was defined as the total horizontal shear force acting between the 
left end of the deckhouse and the section at z. The total shear force 
in the structure V is equal to dM/dz. According to Equation 3, the 
shear T at the ends of the deckhouse is not zero. At points slightly 
away from the deckhouse ends, there is no deckhouse, and hence T = OQ. 
This means that in addition to the distributed shear force dT/dz, there 
must be concentrated shear forces i and Ty at the ends of the deck- 
house. These concentrated shear forces cannot exist in any real struc- 
ture. Their appearance is due to the fact that shear-lag effects were 
neglected when Navier's hypothesis was assumed to be correct for the 
full length of the deckhouse. These forces will, in reality, distri- 


bute themselves over a finite distance, presumably over a distance equal 


to the depth of the deckhouse as depicted in Figure 1.3. 





FIGURES]. 3 








With this in mind, the stresses found using Equations 2 will be 
incorrect for calculations near the ends of the deckhouse. Accord- 
ing to Saint Venant's theorem, however, the effects of the simplifica- 
tion will not affect the stresses in the center portions of the struc- 
Bure’, 

In simplifying the expressions for the moments and direct forces, 


the following notations are used: 


teat Le (5) 

A Ay a7 Ay 

eae he (6) 
I. a a5 I, 


The term j!) is a measure of the size of the deckhouse in relation to 


the hull and is referred to as the size factor. The total moment of 


inertia I of the section can be expressed in terms of Ty. I, and by 
the constant I, 
met. + 1. + I : (7) 


A 2 A 


With the use of notations (6-7), the following important expressions 


are derived: 
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These expressions, along with Equations 2, could be used to determine 


0; and 0. in the deckhouse and hull. It will be noticed, however, 


that Equations 8 contain two terms the first term being the value of 
the respective N or MN if the conventional bending theory were applica- 
ble to the entire section. The actual stresses O can now be expressed 


as the sum of the stresses On according to Navier's hypothesis and a 


correction (Ao). 


Oo = 0, + Ad (9) 


Math the wavicr's stresses (3...) 
a 


a, = - = (10) 


where x is the distance from the centroid of the entire section, 
Figure 1.1. The corrective stresses 40, and Ao. in the deckhouse 


and hull, respectively, are 


Ao. = aL = — x (lla) 
1 AL I, = 
AN, AM, uy 

Ao, = i, = T x (11D) 


where 4:1] and A‘{ are the corrective portions of N and “ given by the 
second terms of Equations 8. 

Although the model in this analysis was simplified in assuming 
Baaeeno vertical forces act between hull and deckhouse (i.e.. zero 
deck stiffness), the important result remains that the stress distribu- 


tion can be expressed as the sum of the stress 0,,, according to Navier's 


N 


theory and the corrective AO. A typical example of stress distribution 


using this analysis is shown in Figure 1.4. 
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In this analysis the assumptivcn is made that any relative displace- 
tnt of the deckhouse with respect to the hull is resisted by the inter- 
tl vertical forces required to deflect bulkheads and transverse beams 
¢pporting the deckhouse. CUxternal vertical loads and bouyancy will 
cuse the structure to deflect. and this deflection can be described 
t the displacements Vy and Y5 of the center lines of the deckhouse and 


ll, respectively, Figure 1.5. 


DEFLECTED CENTERLINE §$ 





C.G. OF DECKHOUSE 








€.G. OF HULL 


FIGERE 1,5 


TM stiffness of bulkheads or deck beams resisting relative vertical 
d: placements is assumed to be constant for the full length of the 
de khouse,and the magnitude of this stiffness is given by a spring 


ccstant K. K is defined by Bleich as being the force per unit length 





of deckhouse required to produce a relative deflection equal to one 
unit of length. The vertical reaction between hull and deckhouse 
will be therefore equal to K(yy - Y>) per unit of length. 

Using the theorem of stationary potential energy, the differen- 
tial equations for the deflections Yy and Y5 can be obtained. The 
potential energy U consists of the internal strain energy and the 


potential energy of the external forces. 
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2/2 
= [Syo Je 72 (12) 
Bewelt be a minimum if the variation 
o= 0. (13) 


Using the calculus of variations, the set of two simultaneous equa- 


Irons are derived. 


Z Ly : 1V : a 

E(I, + ayy + Ky, ~ a, a, ELY5 - Kg a Ps (14a) 
iv , , DV af 

a4, EL} - Ky, + E(I, + aol dy, + Ky» any (14b) 


Py and P5 are the distributed loads acting on the deckhouse and hull, 
respectively. 

The process of the calculus of aeons also furnish the boundary 
conditions required to determine eight arbitrary constants which will 
appear in the general solutions of the differential equations. At 


zZ2=+ 2/2 and z = = &£/2 








Om 


= (15a) 
2 ' § 1? 
E(I, + Cy TY, + Ob %5 EI, an 0 (15b) 
. ' ? a ? ? ? , 5 
O14 5 ELV) + E(I, ct, Ly. = — if (15c) 
2 | Tee 0 
E(I, ~ ay Ty + OL Be, = (iad) 
The second and third boundary conditions are equivalent to ary =~ aa,T 
and If, = M - aa,T as presented in Fquations (1), The fourth boundary 


Memeition indicates that the shear forces at the end ot the deckhouse 
must vanish. 


Seeeeocolution of Differential Fauations for Constant ‘oment 
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In reference to the free=-body diagram of Figure 1.6, this analysis 


neglects the loads Py Po and the shears S. and Sy The only load on 


C 


the model is Mo = “1h = If (i.e... constant moment). The differential 
Equations 14 are then homogeneous, and since the problem considered 


is symmetrical with respect to the origin of the coordinate z axis, the 


general solution contains four arbitrary constants. 


2 : 

a C, ac C42 et C. SinyZ *saminy 2: Cy cosy2z coshyz (16a) 
Zz : 

Yo = C, + C,2° - uC, sinyz sinhyz - uC, cosyz coshyz (16b) 

where 
A 
K 1+y 
aa CF (17) 
4E aot, ap 119 


These solutions in conjunction with the boundary conditions furnish the 


values of the constants, 
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[7 «CET 
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4 ay? (1 + 1) E(a,1, + ua, I.) 


Cc) is not included since it only describes rigid body motion. 


values of 1 and Vy are 
Sim WU cosh tu + cos Ww. Gini ou 
i ean COS UW +t Gin t- cos u 


SoS uu Sinn u = sin wu cosh. u 


Vy Sin u cos u + sinh wu cosh u 
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Introducing the equations describing the deflections a and V5 


the following expressions for the moments and direct forces 


"7 


My = - El, yy 
,? 
iS = = | 
M, EL, ¥4 
EL 
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1 e a ie 2 


(18a) 


(18b) 


(18c) 


The 
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(19b) 
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Ineo 


(21a) 


(21b) 


(21c) 


the following values were found at the amidhsips section, z = 9. 


ee ee 
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A remarkable result is that if a comparison is now made between Fquations 
8 obtained in Section 1 and Equations 22, the only difference is that 

the factor y is now introduced into the second term. Therefore, the 
refined theory used in this section furnishes a similar result as that 
obtained in Section 1 (zero deck stiffness) but modified only by the 


multiplication of the corrective stress Ao by a factor >, 


Co = 0, + ?, NO 3 C3) 


is result is surprisingly simnle in that the deviation of the stress 
distribution from Navier's calculations is determined by a non-dimen- 
eaenal factor 21° The corrective forces Any. AN, and moments AM) AN, 


are computed from the equations given in Section l 
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FIGURE 1.6 


4. Solution of Differential Equations for Equally Distributed Loads 
In this final analysis Bleich considered only the case of equally 
distributed loads Py and P5 acting on the deckhouse and hull, respec- 


tively, where the moments at the ends of the deckhouse are set equal to 





Bero (i.e... M. = M. = 0). Equilibrium requires that the external. shear 
forces Sy = “Se = d1. (p, + Do)s The moment at the center of the sec- 


tion (Figure 1.6) duc to the loading is 


® Se 
"\p 8 a (24) 


The general symmetrical solution to [Equations 14 are 


ee c. + Coz" + Cowcaney2 -Sinimey2 C, cos YZ cosh yz 
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The boundary conditions expressed in [Equations (15) furnish the values 


Gime constants 
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and where 


2 sin u sinh u 

Po ~ usin u cos u + sinh u cosh u (27a) 
a cos u cosh u 

V5 ~ usin u cos u+ sinh u coshu ° (27b) 


The values of the constants are similar to those in the preceding 


section, and the remaining solution for the stress distribution at 





=e 4 


Ze=O is the same, the only difference being that a new deviation 
factor %, replaces °) and - replaces 4. 
= +O aS 28 


Curves of % and %, ase plottwsedeas a fwumetion of u in Figure 1.7. 





$ 





FIGURF 1.7 


meee is a function of KK, it will be noticed that for the case of 


zero deck stiffness (KF = 0) that a and 5 are both eaual to 1.0,and 


hence the results of section one are confirmed 9 = 7. + a7. then % 


Or the dirensiors of the deckhouse became laree and hence U > 2, the 
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Although Bleigh's analvsis dealt primarily with the stress dis- 
Mmewition located at the center of the deckhouse (z = 9), it is the 
intention of the present analysis to also investigate the stress 
Uastrivutions away from the center. <Altnourn Pleich mentions that 
the relationship for stresses along the declk:house will be in the form 
of 

o=0., + @(z) to, (1) 
no expressions for ¢ as a function of z ''ere provided. They can, 
however, be readily formulated. 
fee oment Loading 

For the case of constant moment loading,the expressions for the 
deflections are: 


Z 
y= C, + C42 <i C4 sinyz sinhyz + C, cosyz coshyz (2a) 


1 


Z 
y, = C, + C,z2 = uC, sinyz sinhyz - uC, cosyz coshyz . (2b) 


The second derivatives of Yue ¥ with respect to.z can be introditeed 


into the expressions for the moments Mays My and Nye Sas 
tf ' f 
y 5 eal, Yy i = -8El, * (a) 
LN tf tf 
Ny2-N,=—- Gy, + oy, ) (3b) 


yielding the following results: 


"one 


My = a [@, cosyz coshyz - Yy sinyz sinhyz] 
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MI, ine ny ayt, z nia, T>) (4a) 








, = + lo. cosvz coctere — ). siny2 Stuliyz | 


2 I ail jel 
Z 
MT, pee ae (45) 
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t i al 1 1 
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Mee semermressions are exactly in the form of Cautions 22 of Chapter 1 
except that a nov deviation. factor now exists such that 


Oo =o, + 0, (2) Ao (5) 


44 


where 


9, (2) = [0, cosyz coshyz - , sinyz sinhyz] : (6) 


1 i 


In order to modify the argument of the trigonometric and hyperbolic 


2 e ° a 
functions, it will be remembered that u = i . and hence yz = 2u 


oo 


w2]N 


Now 9, (z) can be expressed in terms of a non-dimensional distance from 


the center of the deckhouse. 


2 Z Z ae 
= 2 AY —-- 1 I ae ae e 
9, (2) e, cos 2u F cosh 2u 77 ¥, sin ou 7 sinh 2u 7 (7) 


Pome the aid of FOR BRAN computer me@thods to facilitiec calculations, 
values of 2, (z) were calculated as a function of u and (z/2). A f£&seely 
of curves for different values of u is presented in Figure 2.1. In 
finding the bending stress distribution away from the centerline, no 
additional calculations need be made other than selecting the appro- 
priate deviation factor from the provided curves and multiplying it by 


the corrective stress Ao. 
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2. Equally Distributed Loads 
For the case of equally distributed loads, the equations for the 
deflections are: 


2 
y, = ©, + C42 + C3 sinyz + C, cosyz coshyz + 





iL 2 3 
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(p, + pidz p 

ee et se 
24 EI (I + wk 

ae tC cs C sinl ee Z. l 
¥9 1 oe WC. sinyz sinhy2 UC, cosyz coshyz 
eae . a 
24 LI aa 


kollowing the same procedure used in the preceding section, the expressions 


for the moments and direct forces are: 
_ MI 


‘ly = B= - [%, cosyz coshyz - v9 Simye sinh | 
pest arn Tie ee a 2° (9a) 
‘ ll a + W) (a1 + HO.) oe iE = 
M, = Bre? 4 [O. cosyz coshyz - Y, sinyz sinhyz] 
ae Ae = oi Rae ee. "2 32 (91 ) 
I, (1 + H) Can ~ LO, I) 2 i 
Lind 
ae, a aT pes aA 5 1 es : 
ao) - oT [5 cosyz coshyz -Y, sinyz Simi 
np_A ee ee - i a ny a (9¢) 
a (1+ u)(,1, + va,I,) aa i j 


At first glance these equations do not appear to be in the standard 
form due to the presence of an additional third term. The first and 


third terms of each expression can be combined, however, such that 
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where Me is the moment duc to the distributed loads located at a dis- 
tance z from the centerline. Equations 9 can then be used to express 
the stress distribution in the form of 

So = of + @» (2) AG Cy) 
where 


9, (2) = 2, cos oe = cosh 2u z - v9 Sain) a - sinh Oi r : (32) 


It must be remembered, however, that the Navier stress term is deter- 
mined through the use of the moment CU) at the location in question 
whereas the corrective stress term (40) emplovs the use of the moment 
at the center of the structure oc) 

As in the previous section, computer calculations vere used to 
find 24 (2) as a function of u and (z/fL). A fSmMily of cutves for differ 
ent values of u is presented in Figure 2.2. 

for the case of constant moment loading. the value of ®, (2) at 
the ends of the deckhouse (2/2 = .5) is equal to 1 for all values of 
u. The reverse is true in the distributed load case. 2, (2) is eaual 
to zero at the ends of the deckhouse for all values of u. This is as 
expected, since the moments at ends of the deckhouse must vanish in 
the distributed load model, and, therefore, there can be no ccrrective 


eeress term at that location. 
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m@emeralized solutions. On tite other hand, there is nothins to indi- 
Meee that tie ayplication of linear superposition is not valid. The 


sirmnlifving method of breaking the free-body diagram cf kigure 1.6 
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Memora constant momeiit part and a ‘tistrivuted lov 
aeeieiate enc slemat forces) + @eem rot wheliwtc eruilibrivm@ nor does 
Memeeeer the actual loadinz, shear, and hendine ro®ent diagrans for 
the model when the ttvo loading conditions are then summed together. 
With the assumption that this reasoning also applies to the respec- 
tive stress distribution solutions, the total stress distribution 


(o_.) will be expressed as the linear sum of the stress due to a 
a 


constant moment (0,,) and the stress due to a distributed loading ft. 


ee 13 
7 3 * OD (13) 
Pere 
CG. = oO. + 9, (2)dc,, (14) 
o go. tO ne OS C5) 


Care, however, must be taken such that the appropriate moments are 
used in these equations. ae for the deckhouse, as an example, employs 
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as described by Equations 11 in Chapter 1. or is different only in 


a¢ 


Mat it makes use of the constant moment ‘“{ in lieu of se CoG oats 


in mind, Equations 14 and 15 may be combined such that: 
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In later chapters this assumption of linearity will be checked by com- 

parison of results on similar models using finite element methods and 


Bleich's method. 
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The Structural Design Language (STRUDL) developed at the M.I.T. 
Civil Engineering Systems Laboratory is a series of computer programs 
for solving structural engineering problems. It is a subsystem of 
the Integrated Civil [ngineerinp System (ICES) and can be applied in 
Many ways to a wide class of practical problems. An operational 
system has been implemented on the IPf system /370. When entered 
into the computer, the STRUDL language is translated by the ICES 
Command Interpreter, which, for each command, calls upon the appro- 
priate STRUDL program. In essence it is a language with which an 
engineer can describe a problem, its solution procedures, and ask for 
results. 

Of the various analyses offered in STRUDL, the stiffness analvsis 
was employed for the present investigation. The stiffness analysis is 
a linear, elastic, static, small displacement analysis. 

The basic principle behind finite element techniques is the 
replacement of the actual physical problem by a model composed of a 
finite number of discrete elements which are connected at their nodal 
points. Many elements are available which have different geometric 
shapes and different applications for the solution of different types 
of meebiens. For the present application, the primary element type 
used is the 'PSR' element. The 'PSR' is a rectangular element with 
four nodal points. one located at each of the corner joints. Used 
in conjunction with a plane stress problem, the nodal unknowns or 


eiepiacements are U, and U,. 
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fomelement loads are available, and hence all loads must be in-plane 
joint loads. (thile resulting displacements are given at the nodal 
points, stresses and strains are given only at the baricenter of each 
element, Figure 3.1. 


fr 





PL CF 


FIGURE 3.1 


The material specified must be isotropic, and the only relevent ele- 
ment constants that need be declared are the Yfoung's modulus and the 
Poisson coefficient. 

Another element type used was the "BPR' element. It is used 
in plate-bending type problems and is similar in geometry to that of 
Seem ese’ element. Tite 'DPR' element allows distributed type loads 
in the out-of-plane direction. In later analyses the 'BPR' elements 
were superimposed on the 'PSR' elements in order to obtain a more 


refined analysis. 
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For this analysis a model vas selected with arbitrary dimensions 
as shown in Figure 4.1. Bulkheads are placed at equally spaced dis- 
tances of 20 feet in the hull section. The thickness of the hull, 
box girder plating is .5 inch, and the thickness of the deckhouse 
plating and bulkheads is .25 inch, The material constants include 
a roung's modulus of 30 x 1 psieeonda -Peaseson cocti fiicrenmyot 0x2. 

ne model is assumed to have a li-foot draft with a corresponc- 
jar mydrostatic unward distributed force of 17.1435 tons/ft. The 
internal loading is arranged so as to provide equilibrium and a 
meeultant symmetrical loading. Figure 4.2. Shear and moment diagrams 
for the total model are also provided. 
1. STRUDL Model 

Due to the symmetry of the loading, it is necessary only to 
look at one-quarter of the deckhouse and annropriate hull structure. 
amis also aids in reducing commuter time which can hecome excessively 
high when dealing with finite element methods. The planes of symmetry 
Seem pronlem arc t've aniéships nlane end the centerline slant of 
the model in Figure 4.1. The one-quarter structure, broken into 122 
See cléements, is shown in Figures 4.3 to 4.5. 

eemee the ‘PSR’ clerents only alloy in-plane joint loads, it sas 
necessary to model the shear, moment, and distributive loads into 
forces acting on joints or nodal points. The pressure loading acting 


@ecie bottom of the hull structure (319.984 psi) was translated into 


EOrces acting at bulkhead locations. The moment at the end of the 
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LOADING SHEAR AND MOMENT DIAGRAM FOR 
MOOEL WITH BULKHEADS 


FIGURE 4.2 
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ONE QUARTER STRUCTURE - ELEMENT ANDO NODAL NUMBERING DIAGRAM 
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decl:house (i) Was tranertated into forces acting on the end joints 
mmmene hull girder. Using the notation in Section 3.3 of reference 


6, the vertical shear flow (°) at the end of the structure was found 


using 
ay, 2 
at of a(l = —) ] 
we (5 2) — —--—___*— = 
2 de 3Af + a 
where 
QO= 5 


D 
meeesmear flow assumed a parabolic distribution, the maximum occurring 
at the neutral axis of the hull. The shear flov was then translated 
into vertical forces acting at nodal points at the end of the hull side. 
Tie final loading is depicted in Figure 4.6. 

A sample computer program for this loading is provided in Appendix A. 
fee = Results 

In order to find the magnitude of any distortion due to the applica- 
tion of point loads on the structure, the variation of longitudinal 
Stress along the deckhouse length was plotted, Figure 4.7. brom these 
results it appears that some distortion occurs near the deckhouse end 
as evidenced by the abrupt change in the stress pattern. These effects, 
however, appear to be concentrated only within a distance of 10 feet 
of the deckhouse end. 

As previously mentioned, the stress output in STRUDL is given at 
the baricenter of each element. It will he noticed, however, that the 


baricenters of the elements adjacent to the amidships’ location lie on 


points 5 feet awav from amidships. There is little or no stress 
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Variation in this area, and hence neplipible error is introduced if 
Sm@ectress at the adjacent elements is used for the anidships' results. 

The longitudinal stress distributions at amicy -ns and at a 
point 35 feet from amidships are shown in Figures 4... and 4.9, res- 
pectively. Pronounced shear lag effects are evident in both cases. 
This is primarily due to the relatively low declkhouse length (1) to 
mele breadth (B) ratio. for this particular model 1/B = 2.5. 

An unexpected result is the discontinuity in the longitudinal 
Stress distribution across the main deck occurring at the junction 
of hull deck and deckhouse side. <An apparent cause for this dis- 
continuity is the nature of the relative deflection of the hull- 
deckhouse connection with respect to the hull side, Figure 4.10 and 
feel, “the longitudinal stress due to bending is proportional to the 


curvature of the beam ) 


g = - E(2) y (2) 
p 

where 0 = radius of curvature. 

It is apparent that the curves of deflection of the main deck 
at side and the hull-deckhouse connection as presented in kigure 4,12 
are quite different from one another. Obviously the curvatures of 
Che two deflections are of different sign or magnitude at respective 
locations and hence effect the stress distribution across the main 
deck so as to produce the results shown in Figures 4.8 and 4.9. In 
Order to support this reasoning. a separate STRUDL program was run 
in which the amidships'plane was held rigid. Although this was an 


unnatural condition, it had the effect of increasing the deck stiffness 
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in the immediate vicinity of the anidshins' plane. The hull and 
deckhouse, therefore, tended to deflect in the same manner.and the 
discontinuity was virtually eliminated as evidenced in tigure 4.13. 

Bmaotnicer area of interest is the distribution of the shear stress 
alone the hull-deckhouse connection. Due to the nature of the ele- 
memes used, it was only possible to obtain results of the shear- 
force distribution 2.5 feet ahove the hulledeckhonie connection. 
Homemer, the trreeular pattern of the shear stress as shoim in kig- 
ure 4.14 is sirilar in forr to results obtained ov exnerirents on 
physical models as descrihed in reference 2. 
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Seemcdagram of bigure 1.6 into a constent mement part and a die- 


Meee) loadine aprt. Wit!. reference to kipure 4.7, 2t cur me 
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Peer ae crCss ctrueture arbvibed 2n }ictire =.1. 2t ta Alrcela tool, 
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ana 2 in Oordaer to arrive at a solution. 

In a model with bulkheads, however, the nrinary difficulty 
is involved with the determination of the deck stiffness or spring 
constant .(K). Since K was defined as the force per unit length of 
deckhouse required to produce a relative deflection equal to one 
ime of length, it is apparent that the value ot K will, in reality. 


vary along the deckhouse length due to the emplacement of structural 


bulkheads. In order to simplify the use of Bleich's method, however, 
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Semverage vValuc of i must be determined. To achieve this and a 
Symmetrical portion of the hull structure was modelled to include 
femeenead and attached deck and bottom plating, tipure 4.15. This 
I-beam type structure is simply supported on its ends and allovred 
to feflect under vertical forces applied to the hull-deckhouse con- 
nections. In the analytical procedure used for determining the 
forces needed to deflect the hull-decknouse connections 1 inch, it 
must be kept in mind that although the deflection due to shear forces 
are negzlipible in most cases, in short deep metal beams the deflec- 
tion caused by shear may become a significant portion of the total 
deflection. In this case shear contributes a major portion to the 
fee deflection. 

The equation for the deflection due to the internal moments (Y,,) 
was calculated through the use of singularity functions. The expression 
for the deflection due to shear ) was obtained through the use of 
Ehe method of unit loads as described in reference 7. The total 
deflection (Y,) was then expressed as the sun of t and ae and set 
Sea to 1 inch. The equation was then solved for P (the distributed 
load acting on the hull-deckhouse connections). The expression for 
K is twice the value of P 

K = 2p, (3) 
This analytic approach to K yielded a value of 1.86 x 10° Dou tor 
the present model under consideration. Sample calculations are 
provided in Appendix B. 

As a further check on K a STRUDL program using 'PSR' elements 


on the model presented in Figure 4.15 was run. Arbitrary forces (f) 





SYMMETRICAL STRUCTURE FOR EVALUATION OF K 


FIGURE 4.15 





were applied at the locations indicated, and the resulting deflec- 
tion at the hull-deckhouse connection noted. The force was then 
Scaled for a deflection of 1 inch. P was obtained by dividing } 
eeeewe width of the flange (240 inches). Using Equation 3,K was 
found to have a value of 1.53 x 104 psi. The disparity betveen the 
STRUDL K and the analytical K can be attributed to the fact that in 
the analytical approach the deflection calculations apply to the 
neutral axis of the team only. 

Sample calculations for the solution of the stress distribu- 
tion at amidships using Bleich's method is provided in Appendix C. 
4. Comparison of “lethods 

Comparison of Bleich and STRUDL results at amidships, 15 feet 
from amidships and 35 feet from amidships, are shown in kigures 4.16, 
4.17, and 4.18, respectively. The results for the stress distribu- 
tion in the vicinity of the deckhouse center compares more favorably 
meemecdo the results at a location 35 feet “yay. This is duc to the 
fact that Pleich's solution neglects shear-lae cffects, as discussed 
in Section 1 of Chapter 1. At each location the Sleich stress dis- 
tribution across the main deck and hull bottom approximates an average 
value of the STRUDL distribution. The Bleich distributions along the 
hull side produce lower peak values of stress as compared to STRUDL, 
The distributions of the stress along the deckhouse side and top 
attributed to Bleich's method are, in general, of greater magnitude 
than that of the STRUDL results. The stress distributions for both 


methods were transformed into moments and checked against equilibrium 
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Semditions: that is, the values of the moments obtained at each sec- 
tion corresponded to the respective moment on the bending moment 
diagram. 

Because of the large shear-lag cffects found in this particular 
model, it is difficult to make an accurate staterent as to the 
applicability of linear superposition to Bleich's methods (which 
ignores shear-lag). It appears, however, that if these shear-lag 
erfects were reduced, equilibrium would reouire that the STPUDL resuits 
approach Bleich's solution. This will be later checked through the 
use of a model whose increased 1/B ratio sienificantly reduces the 


shear-lag. 


CHARP ey =< Ope ITT MCh GASTD 1/8 


In this analysis the same model as presented in Chapter 4 was 

used except that all dimensions in the longitudinal direction were 
M@ienceascd by a factor of to so as to increase the 1/3 ratio from 
memto 2.9. %In order to maintain a bendine-moncnt diagran similar 
to the one used in the preceding chanter, the distributive loading 
was reduced hy a factor of four: i.c.. p, + Py = 1.4285 Coneyert . 
From this lording So amo = Loa, ance = “eo a = / Gs? @ pevcons. 
Due to the increased distance between bullheads, the valuc of the 
Meck stiffness (!:) was refuced to 7.67 x 197 pean. 

The sare procedures as described in Chapter 4 vere used in the 


determination of the STPUDL and Bleich solutions. 


1. STRUDL Results 





The longitudinal stress distributions at amidships and at a point 
mmereet £Lrom amidshins are shown in Figures 5.1 and 5.2, respectively. 
The most netible dcifference in these results from that cf the precec- 
ing chapter is that the deckhouse has become more effective in the 
Suoport of the stress distribution. This is due to the fact that the 
deckhouse leneth has been increased to 209 feet. 

The discontinuity in the stress distribution across the main 
deck is still apparent but to a lesser degree. <A look at the deflec- 
tion curves of Figure 5.3 indicates that the differences in curva~ 
tures between the deckhouse connection and main deck at side are 


not as great as in the preceding model. 
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fiyeemast imeortant result. however, is that tt’ shear-lac effects 
have been reduced to a neglipible anount. 
2. Comparison of tethods 

Comparison of STRUDL and Bleich results at amidships. 30 feet 
from Madiships and 79 feet from amidships, are shown in Figure 5.4, 
5.5, and 5.6, respectively. Again, the difference bettveen solutions 
grows progressively worse as one travels away from amidships or the 
center of the deckhouse structure. 

The most imnortant result of this analysis is that at the amtd- 
ships section (Grhere Bleich's solution is most applicable). the two 
solutions are almost identical. The argument that linear superposi- 
tion is applicable to Bleich's two general solutions is therefore 
confirmed by these results. It is readily evident that if the shear- 
lag effects were completely done away with, the tvo solutions would 
correlate exactly, except for the discontinuity at the hull-deckhouse 
Memeetdo:, Li must be képt in mind, however. that this savs notfing 
Or the applicability of Bleich's total solution for use as a design 
tool, In real life structures shear-lag effects are inevitably 
present to some degree as exemplified by all the STRUNL results thus 
far presented. Further discussions on this topic are presented in 


later sections after all model results have been presented. 
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CHARGER VI = MODEL VITEOUT BULIGIZADS 


In order to investigate the effects on the stress distribution 
due to a more flexible deck, the bulkheads of the model presented 
in Chapter 4 were climinated. All other parameters such as loading, 
lengths, boundary conditions, etc. remained unchanged. 

kor the STRUDEL model, '@SR' and 'BPP' elements were super- 
imposed upon the tole structure so as to allow out-of-plane loads. 
Mento the nature of the iwaposed boundary conditions, which allows 
a considerable degree of flexibility, the program was not acceptable. 
It is also noted here that for models with bulkheads and using both 
"PSR' and 'BPR' elements, the costs were prohibitive due to the 
increase in the number of elements. In order to obtain reliable 
results for the present analysis, only 'PSP' elements were used. The 
loads on the hull bottom were attached to transverse stiffeners in 
the same locations where the bulkheads had once been. 
i) STRUDL Results 

The longitudinal stress distributions at amidships and at a 
noint 35 feet from amidshins are sho-m in Figures 4.1 and 46,2, 
respectively. The most notable result is in the degree of stress 
reversal occurring along the hull and deckhouse sides. Because of 
tne relatively flexible deck, the deckhouse has hecore less effec- 
Give in the support of the stress distribution. 

Shear-lag effects are as pronounced as in the original model: 
however, the discontinuities at the hull-deckhouse connections are 


not as great. Figure 6.3 supports the diminution of discontinuity 
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as well as hiphlightine the increase in relative vertical dienlace- 
ments between the hull-deckhouse connection and main deck at side. 

inee sheAr-strass distributiow 2.5 feet mhove the hull-deck- 
femee Connection is showm in Figure 6.4. 


2.) Baeich Model 


ee 2.9 ane ~ * ae ee 


The absence of bulkheads facilitates the analytical approach 
Seeeie da@termminetion of K. A transverse slice of the main dec! 
@aarebe assumed to he a simpnle slender tearm rigidly svu-~nerted on ite 


ope «6 lfc hear ie erhiect to vertical deflections due to distrthuted 


LN 


Joads (P) so as to deflect the hull declhouse cornnecticrn 1 ine’. 


eam, simple-beum wicilection theory. the valuc of K was found to 


Pel. cwididevelii is very sesacl wiyen cowmamed cv tle “ for 


é 


ne love 
m@ilotner bulichead models. 


3. Comparicon Sf Met':ocs 


ate 
OR A ep a EE 7 em = ae ee eet 


TY V IIS 


Gemparison of Bleich and STUURE results at amidships. 15 feet 


Poe mee Cellar 27d @&5 feet fram, vwet dehi-k. 2re cheer in |] ifuree 6.5 


eee 0./, respectively. MRin. the diffard@nee Betrecn solutions 


Peo eseproecoSssive.y worse as one travdls awey fr @m tic Weidgeins area. 
ane Bleich dist:ri>uticns along the null side produce lower peak values 
of stress as comnared to STRUPDL. whereas the distributions of stress 
along the deckhouse top and side are, in general, of greater magni- 
tude than that of the STRUDL results. In each case the Bleich stress 


distribution across the main deck and hull bottom approximates an 


average value of the STRUDL distribution. 





S SEN 4 


t 
& 
t 
=a 
’ 
. : é 


~ OlX ISd Nt SS3NLS 


4 
( 


NOISNIL | NOISS3YdWO) © 
Of SP? OF G) oO} S ° cS 0} Ci or Ge oO£f SE€ 





~ 


— 





i 

; | ; 
a 
1 TSAD7 N70 NIVI - i 

| 
-_—_ -—- | ee ee 
VWONYIS 
941 =: } 
L SAIHSQIWy LY "Ga SS Fees FS 


Gt =8/¥ SavaHrNG Sr Lun “13.00 | 
| , 


> SGOH1aW SO“NOSI¥VdWOD. ~ 


aoe ee. ined yi | ae 








i 


seer] Ft 2 OlX ISd 





@eo+e 


Sal HSOlWW | Wold 1334 S$) 


=i 
‘ 
' 
; 


, 9 OVUMNDL 4 


~ 


$d 


tp de 


Ge =9/%, $,0HE LNOULI 


eh | 
W 4JO.NOSI4 


| 
O 


S€ Of Ge or Gi Ol ¢ 


V3A31YI30 NIVW 


413 


rrr 


NOISNAL 


em 


El 4S ae 
YGANW1S 9] ae 3 
1610 SSaNIS 
T30d0W a : | 


| 
| 
oO ¢ 


3 
Vd 


ey 


mM - 


| 
Rt ssauls 


| 


| 


es 


| | 
bl 
i fer | t ai i | 2“ 





! 


NOISS3IYd 


WO?) 





Oe St Of GE 


| 
| 


iz 








£°9) FNL A | | 


~IIXISd Ni SSavis 





NOISN3L | NOISSAYdWOD 
Of S$ CX a5 Bot os Oo S o18 St for srijo¢ 


| 
| ‘ : 








| 
HDI1Ee— —-— ~CO 4 
| VGNYLS | | 
7 oe el) el o ee ICL = Hy 
“SAIHSQIWNY WOUs LAR4 SE “1910 SSAHLS | 
ea | | | — Se = Sh SQHE LOOHIIM - -1300W 


aes 


te} ES) |i caguianao FNosiuwaWwO2 


a at ete al = 


her —— nen — it — ae eee 





=i 
CORCLUS tS ARM @LOQERE DAT IONS 


The principle of linear superposition can be applied to “leich's 
two general solutions in order to obtain the total stress distribu- 
tion solution. Because Bleich's method neglects shear-lag. the 
use of this method may be of limited significance when dealing with 
Pee iafe structures. As the 1/8 ratio increases. however, Bleich's 
solution becomes more applicable for the center portions of the 
decehnouse—-hull structure. 

pecause ot the relatively favorable comparison of Bleich and 
STRUDL results in the center portions of the deckhouse-hull area for 
all models considered, there is some indication that Bleich methods 
may have some application for design pruposes. The rather simple 
analytical metnods which are employed for the determination of a 
Steess sOlution is also a factor in consideration. All of the results 
of the present analysis indicate the Bleich method is conservative in 
reference to the stress distribution over the deckhouse proper and 
that it supplies an approximate average stress across the main deck. 
The only shortcoming is that the Pleich solution underestimates the 
value of the stress at the junctions of the hull side with the main 
deck and hull bottom. kor example, in the model with bulkheads and 
Wiper 2.), the value of the stress at the main deck at side is almost 
20 per cent less than the STPFUDL result. As the deckkouse length is 
increased so as to diminish the effects of the shear-lag, Bleich's 
solution converges towards the finite element solution, as evidenced 


in kigure 5.4. This supports the conclusion that Bleich's method is 





Boe 


more strongly applicable to structures having relatively long deck- 
houses. 

In order to confirm the finite elerent results. esnecially in 
the area of the discontinuity occurring at the hull-deckhouse connec- 
tion, it is recommended that a physical model be built similar to 
the mathematical models presented in this analysis. Strain and 
deflection analyses under similar loading conditions could then he 
used for verification, and augmentation of results obtained throuph 


SERUDL. 
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The following is a STRUDL program for finding element stresses 


and nodal displacements for a model with bulkheads spaced 20 feet 


apart. 
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The following is an analytical approach to the determination 


of the spring constant (K) for a model with bulkheads spaced 20 fect 
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NOPRRNTIM Com SA PIT BIRICH CALCURATIONS 


mee followime 16 a solution for the stress distribution Wt 


amidshins for a model with bulkheads spaced 20 feet anarct. 
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